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We review some of the features of the evolutions equations for transverse momentum dependent 
parton distributions recently proposed by us. We briefly describe the new ingredients entering 
the equations and their relationship with ordinary evolution equations. We comment on possible 
choices for the initial conditions and then show results for the evolved distributions obtained by 
numerical implementation of the equations. By computing the average transverse momentum at 
different scale we highlight some general properties of the solutions. 
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1. Introduction 

Transverse momentum dependent, or equivalently, -distributions are currently object of in- 
tense research activity due to their wide range of applicability in the description of hadron initiated 
hard processes [|l], The use of /c^ -distributions is indeed phenomenologically appealing since 
observables constructed upon them show a reasonable agreement with data already in lowest or- 
der, which is not the case for predictions based on coUinear factorization at the same accuracy, see 
for instance the discussion in Ref. Their correct formalization in quantum-chromodynamics 
appears to be well established only at high energy Away from this limit the situation is still 
unclear, although detailed investigations are present in the literature [Q, ^, The relevant issue 
of factorization within a -dependent approach has also been investigated [^]. The latter prop- 
erty together with a precise knowledge of the scale dependence of -distributions would allow to 
relate to each other data coming from experiments at different energies and, more important, to 
test factorization quantitatively. Although a definitive answer is absent in the literature, there have 
been however some attempts, see for example [I^. In particular we will focus on the equations 
proposed by us in Ref. JTl] ] which are the space-like version of the ones proposed in Ref. [|l^]. In 
a subsequent phemonelogical study [[T3|], performed in the context of semi-inclusive deep inelastic 
scattering, it was shown that a reasonable description of data could be obtained once unintegrated 
evolution equations were solved with suitable, but motivated, initial conditions and assuming fac- 
torization for the cross-sections of interest. This result stimulated us to apply the same formalism to 



the description Drell-Yan type processes in hadronic collisions [14]. The -spectrum of the gauge 
boson has, in fact, a rich transverse structure and involves both perturbative and non-perturbative 
aspects of the underlying theory. The resummation of the perturbative series in the multiple soft 
gluon emission limit can be accomodated by using properly modified unintegrated evolution equa- 



tions [15] so that the structure of the non-perturbative form factor can be investigated. 



2. Space-like fcf -dependent evolution equations 



Transverse momentum dependent parton distribution function, ^'p{xB, Q'^,k±), give the prob- 
ability to find, at a given scale Q^, a parton / with longitudinal momentum fraction xb and trans- 
verse momentum k± relative to the parent hadron momentum. The evolution equations for fc,- 



distributions proposed in Refs. [11, 13 1 read 



Q 



2k 



du 



dH^ 
n 



Xb 



(2.1) 



and resum large logarithms associated to the emission of coUinear partons. However, at variance 
with ordinary evolutions equations [|T6|], the transverse momentum generated at each branching, l^, 
is explicitely taken into accont kinematically and eventually adds to the non-perturbative one due to 
Fermi motion of the parton in the parent hadron. If we now consider to be the bare distribution 
of a quark in a quark, neglecting flavour indeces, we may insert the source term ^(x, = 
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Figure 1: Pictorial representation of the proposed evolution equation. The longitudinal and transverse 
component of the partons involved in the splitting k ^ k + l ax& indicated. 



5(1 — x)5^'^\ki_) into eq. (2J_) obtaining, up to ^{a^), the following result 



^(x,G2,*^) = 5(l-x)5P)(^^) + gl^p(x) 



(2.2) 



which indeed exposes the singularity associated to the emission of a collinear parton weighted 



by splitting functions P{x) [p^. In order to clarify the role of the new ingredients appearing in 
eq. (2.1), it proves useful to consider just one emission, see Fig. (|]). The small blob at the bottom 
left of the figure then symbolizes the iteration of emissions in the parton ladder of whose only 
the last, ^ — )• ^ + /, is explicitely shown. Transverse momenta k± and k± are defined with respect 
to the incoming hadron direction, while l± is defined with respect to k direction. The transverse 
momentum k± and longitudinal momentum x of the branching parton k are written as a function 
of the relative transverse momentum l± and the fractional momentum u appearing in the splitting. 
They can be expressed as k± = {k± — l±)/u and x = x/u, respectively. The former is a result of 



the transverse boost from k to hadron direction []17p. Both can be found in the distribution under 
convolution in the right hand side of eq. (2.1). The additional integration (i^/^and the 5-function 



appearing in eq. ( [2.1| ) are associated with the phase space and mass-shell contraint of the emitted 
partons. The latter can be written as 



I' 



[\ - u)k^ + u{\ - u)k^ , 



(2.3) 



where the emission of massless partons has been assumed (/^ = 0). As is well known, logarithmi- 
cally enanched contributions arise when the interacting parton virtuality increases along the ladder. 
Therefore in the hmit k^ ^ k^ we may neglect k^ with respect to k^ and set —k^ = Q^, obtaining 
/2 = (1-M)g2 which can be found in eq. (^). Furthermore, we assume ^, -distributions to fulfil 
the normalization condition: 



d'k^^l{xB,Q\k^_)=f'p{xB,Q' 



(2.4) 



where fp are ordinary parton distributions. It is important to remark that integrating eq. (2.1) over 
(f-k±^ and by using eq. (2.4) we recover the ordinary evolution equations for /p [16] . 
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Figure 2: Left and middle panel: average transverse momentum (fe^) for up-quark and gluon at three 

different scales: Ql = 5 GeV^ (-), = 10 GeV^ ( ) and = 20 GeV^ ( ). Right panel: the 

transverse spectrum of the up quark distributions at fixed = 0. 1 at the previously indicated three different 
scales. 

3. Numerical solutions 



The unintegrated evolution equations, eq. (gjj), are numerically solved by a finite difference 
method on a discrete grid in {xb , ) space. The initial conditions at the starting scale Q^ = 5 GeV^ 
are chosen to have the factorized form 



1 717 



^l>{^B,Qiki_) = fp{xB,Qt;) 72 e^'^^' i = q,q,g (3.1) 



— e 

l2\ 



where /p(jcb,2o) ordinary parton distributions [18]. The choise of a energy (xb) independent 
guassian transverse factor in eq. ( p| ) is supported, for example, by recent semi-inclusive deep in- 
elastic scattering data at low energy [|T9|]. The width {k]_ .) are fixed to a value of 0.25 GeV^ both 



for quarks and gluons ||13|]. Since the initial conditions, eq. ([3.1|), by construction fulfil eq. ([2.4|), 
the latter can be checked for every values of xb and to estimate and, eventually, to improve the 
numerical accuracy of the evolution. For this check to be meaningful we adopt the same flavour 
scheme and coupling values used in Ref. [|T8|]. The evolved distributions show indeed some in- 
teresting properties. In the left and middle panel of Fig.@ we present the average transverse 
momentum, (^^), as a function of xb calculated at three different scale for up quark and gluon, 
respectively. In both cases, the averaged transverse momentum at the final scale, = 20 GeV^, 
increases with decreasing xb- The effect is more evident in the gluon case, due to the singular 
behaviour of Pgg{u) at small u. In the opposite limit, xb ^ I, the averaged transverse momentum 
approaches the value given in the initial condition and indicated by the horizontal line in left and 
middle panel of Fig.(^. In this limit, phase space only allows the emission of soft partons which 
generate negligible trasverse momentum. In the right panel of Fig. (|2|) it is shown the up-quark 
distribution at fixed .xg = 0. 1 as a function of k\. Interestingly, the factorized form of eq. ( p7T| ) 
is not preserved under evolution. The guassian dependence on k\ at is turned into a inverse 
power law 1 /{kY)" at the final scale Q^, the latter dependence in qualitative agreement with matrix 
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elements behaviour at large k^, eq. (2.2). We also observe that, as Q increases, the distributions 



start to populate the higher k\ region, as an effect of the leading logarithmic approximation built-in 



the evolution equations. 
Conclusions 

In this contribution we have reviewed some of the features of the unintegrated evolution equa- 
tions proposed by us and discussed some peculiar properties of the evolved -distributions obtained 
by numerical solution. The use of unintegrated distributions is indeed phenomenologically appeal- 



ing [^, |l4p. Most interestingly, the analysis of jet observables within this framework may give 
access to valuable informations on the unintegrated gluon distributions which describes hadron 
structure as probed in high energy collisions. However, the departure from a pure coUinear scheme 
introduces a number of non trivial problems which are at present under investigation by many 
groups. This implies that, to date, any conclusive statement on the validity of ^, -distributions ap- 
proach in general, and on the proposed evolution equations in particular, must wait for further 
theoretical developments. 
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